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Abstract 

Using a statistical model for the normally deformed states and for their coupling to a 

member of the superdeformed band, we calculate the ensemble average and the fluctuations of 

the intensity for decay out of the superdeformed band and of the intraband decay intensity. 

We show that both intensities depend on two dimensionless variables: The ratio T^/Ts and 

the ratio T^/d. Here, I> is the spreading width for the mixing of the superdeformed and the 
O ! 
Q\ . normally deformed states, d is the mean level spacing of the latter, and r$ (Tn) is the width for 

Q\ ' 

gamma decay of the superdeformed state (of the normally deformed states, respectively). This 

parametric dependence differs from the one predicted by the approach of Vigezzi et al. where 

the relevant dimensionless variables are T^/Ts and Fyd. We give analytical and numerical 

results for the decay intensities as functions of the dimensionless variables, including an estimate 

of the error incurred by performing the ensemble average, and we present fit formulas useful 

for the analysis of experimental data. We compare our results with the approach of Vigezzi et 

al. and establish the conditions under which this approach constitutes a valid approximation. 
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1 Introduction 

The study of super deformed (SD) bands is one of the most active fields of nuclear structure 
studies at high spin |I| . The intensities of the E2 gamma transitions within a SD band show a 
remarkable feature: The intraband E2 transitions follow the band down with practically con- 
stant intensity. At some point, the transition intensity starts to drop sharply. This phenomenon 
is referred to as the decay out of a SD band @ . It is attributed to a mixing of the SD states 
and the normally deformed (ND) states with equal spin. The barrier separating the first and 
second minima of the deformation potential depends on and decreases with decreasing spin /. 
Decay out of the SD band sets in at a spin value Iq for which penetration through the barrier 
is competitive with the E2 decay within the SD band, see Ref. ||. 

The theoretical description of the mixing between SD and ND states uses a statistical model 
for the ND states first proposed by Vigezzi et al. j|, [| |5| . The ND states to which the SD state 
is coupled, lie several MeV above the ground state. The spectrum of these states is expected 
to be rather complex. It is, therefore, assumed that the ND states can be described in terms 
of random-matrix theory or, more precisely, by the Gaussian Orthogonal Ensemble (GOE) of 
random matrices J3J. Likewise, the El decay of the ND states is calculated within the statistical 
model. The results of this approach have been used to analyze experimental data |7|, ||, ||, [10]. 



It is one of the aims of such work to determine the strength of the coupling between the SD 
and the ND states and, thereby, properties of the barrier separating the first and the second 
minimum of the deformation potential. We mention in passing that recently, experimental 
evidence for non-statistical decay out of the SD band has been found |TT| in 60 Zn, a nucleus 
very different from the nuclei investigated earlier |7|, || |9|, |1| . 



The model of Ref. || was re-examined in Ref. |[L2|| . This was done because inspection shows 
that the formulae derived in Ref. || apply only in the limit where the electromagnetic decay 
widths Tjy and T$ for the ND and the SD states are small in comparison with the spreading 
width T^ which describes the mixing of the ND and the SD states. However, analysis of the 
data |7|, ||, || |l(J has shown that this condition is not met in practice. In the present paper, 



we follow up the work of Ref. []12] in which attention was focused on the ensemble average of 



the intraband decay amplitude. We evaluate the contribution of the fluctuating part and show 
that it cannot be neglected in situations of practical interest, in contrast to the claims made in 



Ref. JT2[. We use the supersymmetry approach developed in Ref. |T3|- With d the mean level 
spacing of the ND states, we show that the intensities for intraband decay and for decay out 
of the SD band depend on the dimensionless variables r^r/r^ and Tyd. We give analytical 
and numerical results for this dependence as well as fit formulas to facilitate the analysis of 



experimental data. We compare our results with those of Refs. [|j] and JT2 



2 Model 

We denote the first SD state with significant coupling to the ND states during the E2 decay 
down the SD band by |0); its energy by Eq; the ND states having the same spin as the state |0) 
by \j) with j = 1, . . . , K and i{>l; their energies by Ej. The ND states decay by statistical 
El emission. We assume that the total El decay widths of all ND states have the common 
value Tn. The matrix elements Vqj connecting the SD and the ND states are responsible for 
decay out of the SD band. This situation is illustrated in Fig.l. We assume that the ND states 
can be modeled as eigenstates of the GOE. Then, the energies Ej follow the GOE distribution, 
and the Vo/s are uncorrelated Gaussian distributed random variables with mean value zero and 
common variance v 2 . The spreading width 1^ is defined as 1^ = 2irv 2 /d. The limit K — > oo is 
taken at the end of the calculation. 

The Hamiltonian H of the system is a matrix of dimension K + 1 and has the form (j, I = 
1,...,K) 

H=( E ° V " ). (1) 

\Voi 5 3l E 3 ) 

To H must be added the diagonal width matrix T,sn given by 

*»=- m ( r ; j J ■ <2) 

The effective Hamiltonian H c g is given by 

H e s — H + T, S n ■ (3) 

The intraband decay amplitude has the form 

A 00 (E) = 1S ([E - FeffJ-^ooTs • (4) 

With 7I = r,s, this quantity describes the feeding of the SD state from the SD state with the 
next-higher spin value, and its subsequent decay into the SD state with the next-lower spin 
value. For simplicity, we assume that the amplitudes for feeding and decay are both given by 
75. Similarly, the amplitudes for decay out of the SD band are given by 

A Qj (E) = ls ([E - H cS ]-\ jlN , j > , (5) 
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Figure 1: Schematic illustration of the decay out of a SD band due to admixtures of ND states. 
The level spacings of ND states are exaggerated in the figure. 



where 7^ = Tat. The total intraband decay intensity has the form 

1 



Aw{E)\ 2 dE, (6) 



27iTs j-00 
and the total decay intensity out of the SD band is 

/ out = — -/ Y.\ A ^E)?dE. (7) 

The identity 

/in + /out = 1 (8) 

follows from unitarity and completeness. Except for Sections |] and ||, we will, therefore, focus 
attention on I in . 

Both Ji n and 7 ou t vary with the realization of the ensemble of random matrices. We are going 
to calculate the ensemble average of both quantities, denoted by a bar. This average involves 
an average over both, the distribution of matrix elements Vqj and the distribution of eigenvalues 
Ej. In any given nucleus, we deal with fixed values of the Vqj, and with fixed positions of the 
ND states \j). In other words, any given nucleus corresponds to a single realization of our 
random-matrix ensemble. The question is: How close to the actual behavior of the system will 
the ensemble average be? To answer this question, we also estimate the probability distribution 
of I[ n . This allows us to determine the error incurred by using the ensemble average. 

While it is possible to calculate J in analytically, the calculation of the probability distribution 
of Jin is beyond the scope of the supersymmetry technique. Therefore, we use a two-pronged 
approach, employing both the supersymmetry technique and numerical simulation. We use 
analytical results for I[ n as a test for the accuracy of our numerical simulation which is then 
used to estimate the probability distribution of Ji n . 

3 Example: Perturbative Approach 

Various discussions have shown us that application of the GOE to the decay out of a superde- 
formed band involves conceptual difficulties. This fact has motivated us to include the present 
Section in our paper. In this Section, we present a simplified version of the GOE approach 
which can largely be dealt with analytically, and which is quite transparent. It is based upon 
a perturbative treatment of the mixing matrix elements Vqj. From the outset, we emphasize 
that this perturbative treatment is not justified in the cases of practical interest, and that our 
work described in later Sections of this paper is not based upon such a perturbative approach. 
Thus, the present Section serves a pedagogical purpose only. We wish to exhibit the problems 



encountered when one applies random-matrix theory to a limited data set of a single physical 
system, and the answers one can give. 

We expand the amplitude A j in powers of Vqj, keep the first non-vanishing term, 

A 0j = ls {E -E + ziy2)-%(E - E j + iT N /2)~ 1 lN + • • • , (9) 

and focus attention on the partial width amplitude 

T ,- = V 0j (E - Ej + iT N /2)- l lN • (10) 

All random variables of the GOE reside in jj. Moreover, decay out of the SD band is (aside 
from trivial common factors) governed by the quantity J2j \lj\ 2 - Therefore, we calculate the 
first and second moments of J2j \lj\ 2 as GOE averages. Using the statistical independence of 
Vqj and Ej and performing first the average over the Vqj, we find 



S>,|W£ ? (11) 



N 

^Te~^~fl' 

3 3 

To perform the GOE average over the energies E~, we rewrite this expression identically as 



YhA 2 = v 2 fdE'- ^ ysfE'-EA. (12) 

Y J (E-E') 2 + T 2 N /4 Y 



By definition, J2j <5(£' / — Ej) = 1/d where d is the mean level spacing. Using F^ = 2irv 2 /d, we 
obtain 



EIt.I^^. (is) 

3 

This implies that to lowest order in the Vo/s, we have 

rVr s . (14) 



'out 



The result (|TJ) is remarkable on several counts. First, we observe that J out depends on the 
dimensionless variable r^/r^ and not, as Ref. Q would have us expect, on Tyd and on TAr/rs. 



Second, we find that J out is independent of T^. This fact is counter-intuitive because for any 
given realization of the GOE, the quantity I ont will surely depend on the magnitude of T^. The 
independence is caused by our use of first -order perturbation theory and by averaging over the 
eigenvalues Ej. This average implies the third curious feature of the result (]Tjj) : We smear out 
the positions of the Ej's completely while in the actual experiment, the decay out of a SD band 
will surely depend strongly on where the ND states closest to the SD state are located. This 
expectation can only be reconciled with the result fll4]) if the distribution of J out about its mean 

6 



value ([L4]) is rather broad. If correct, this statement would imply that a determination of T^ 
from experimental data using the statistical approach would necessarily involve large errors. 

To check this point we calculate the variance of J2j \lj\ 2 - A straightforward calculation 
shows that the variance is proportional to (T^) 2 (d/T^). For the variance of J out , this implies 
(-^out — -^out) 2 ~ (T^/T s ) 2 (d/T N ). This result confirms our expectations. For T N <C d, i.e., for 
sharp ND states, the variance of J out is bigger than (J out ) by the huge factor d/T^. The wider 
the ND states are, the smaller this factor. For T^ = d, the factor is unity and decreases further 
with growing T N /d. This is physically plausible: For T^r <C d, the decay out of the SD band 
will sensitively depend on the precise positions and decay properties of the ND states, leading 
to a large uncertainty in the prediction afforded by the ensemble average and, hence, in the 
value of T^ deduced from the data. For obvious reasons this dependence of the decay out of the 
superdeformed band on individual properties of the ND states decreases with increasing width 
of the ND states. 

Clearly, lowest-order perturbation theory in Vqj is not an adequate way of dealing with our 
problem. Nevertheless, we learn from the example of this Section that calculating only the 
ensemble average of I m or of 7 ou t is not sufficient. We recall that in the mass region A ~ 190, 
d/Y^ is typically of the order 10 2 . In this case, we need information on the entire probability 
distribution of I in if we wish to assign an error to results obtained from comparing experimental 
data with the GOE average. This point was also emphasized in Ref. [|J. 

4 Ensemble Average 

Following standard procedure in the statistical approach, we write A 00 (E) as the sum of the 
average part A 00 (E) and the fluctuating part Aqq c (E), 



A 00 (E) = A Q0 (E) + A^ C (E) . (15) 

Calculation of the average part is straightforward JT2I and yields 



A m {E) = ls {E -E + iTs/2 + 1^/2)^3 ■ (16) 

The ensemble average modifies the propagator through the SD state by the addition of an 
imaginary term iT^ /2, This is well known from the theory of the optical model in compound- 
nucleus scattering. The decomposition ( |T5j ) entails a corresponding decomposition of I m) 

T~ = 7^ av + 7~ fluc (17) 

-Mil -Mil T -i m i \ ' ) 



the two terms on the right-hand-side being defined in terms of |Aoo(-E')| an d of I^toC-^)! 



respectively. We have 12 
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We observe that for T^ <C T s , this result agrees with the perturbative result for I out in Section |3|. 
We turn to the calculation of I m . In Ref. |Oj, it was argued that for Tn 3> Ts and 



r^v ^> r^, this term is negligibly small because the ND states decay overwhelmingly by El 
emission. While this assertion is certainly qualitatively correct, the question remains: How big 
is the term I in in comparison with I m , and how does it depend on the parameters Ts, Tn, T^ 



and d of our model? To answer these questions, we use the supersymmetry formalism [1^]. We 



do not reproduce here the complete calculation which is lengthy but quite straightforward. It 



runs parallel to that of Ref. |T3]|. Rather, we describe a shortcut which suggests the form of 
the final result and which also lends plausibility to this final result. 

The formalism of Ref. [jn| is taylored to compound-nucleus scattering. We use the fact that 



formally, the present problem has much in common with compound-nucleus scattering: The 
ND states are compound-nucleus resonances which may decay either by El or by E2 gamma 
emission. The main difference to standard compound-nucleus scattering lies in the fact that 
the SD state |0) acts as a doorway state. Population of the ND states from the SD band, and 
decay of the ND states back into this band, both proceed only via intermediate population of 
the SD state |0). 

To display this similarity, we define the quantity 

Soo(E) = I - iAo Q (E) . (19) 

We claim that Sqq(E) can be viewed as a bona fide S-matrix element. To make this claim 
plausible, we consider first the case where Vqj = 0, for all j. Then, Sqo(E) has magnitude one 
and the form 

This is a one-dimensional unitary S'-matrix describing elastic scattering with a resonance lo- 
cated at Eq of width Ts- For Vqj ^ 0, the coupling of the SD state to the ND states and the 
ability of the latter to undergo El decay, open additional decay channels. Then, the magnitude 
of Soo is smaller than unity, and Soo may be viewed as one element of a unitary S'-matrix 
comprising the El decay channels in addition to the E2 SD band, and displaying the (K + 1) 
resonances stemming from the SD state |0) and the K ND states \j) with j = 1, . . . , K. The 
actual construction of the other elements of this S-matrix is not needed, of course, since we 



are only interested in Sqq(E). In analogy to Eq. (p~5|) we write 

S 00 (E)=^ jE) + S^ c (E). (21) 

From Eq. fllSD we have 

# - E Q - iTs/2 + irV2 



5oo(jE;) " £ - E + l T s /2 + ffn/2 ' (22) 

The fluctuating parts of v4 o(-E) and of Sqq(E) differ only by the factor (— i). Therefore, we 



can use Eq. (8.10) of Ref. @ to calculate |A(JS C (£)| 2 by calculating \S^ 1C (E)\ 2 . The input 
parameters of this equation are specified as follows. The channel indices a, b, c, d are all equal 
to 0. The transmission coefficient T which couples to the SD band is given by 



T = 1 - \S 00 (E)\ 2 



(23) 



(e - E y + (t s + ri) 2 /4 ■ 

This coefficient displays a resonance at E — E with width Fg + TK This is due to the fact 
that the SD state |0) is a doorway state for formation of the ND states from and for their 



decay into the SD band. The parameter e in Eq. (8.10) of Ref. [H| is given by the difference of 
the energy arguments of two scattering amplitudes. The energy arguments of A^i^E) and of 
Aqq c (E)* coincide, suggesting that we put e = 0. However, since the El decay of the ND states 
is summarily accounted for in terms of their common width Tn, an imaginary energy difference 
arises which amounts to the replacement e — > — iT^. As a result, El decay is accounted for by 
the appearance of an exponential factor exp(— (7iT N /d)(\i + A 2 + 2A)) under the integral. All 
transmission coefficients except T must be put equal to zero. The resulting equation expresses 
|Aqo c (-E)| 2 as a threefold integral over real variables A, Ai, A 2 . For the calculation of I[ n , we 
need to integrate in addition over energy E. This yields eventually 

l^ fluc l f°° AT? f°°A\ f°°A\ ^A\ 

Im = , r = / dE / dAi / dA 2 / dA 
lD7rl s J-oo Jo Jo Jo 

(1_A)A|A 1 -A 2 | 

X ■ 



2 



((1 + Ai)Ai(l + A 2 )A 2 )V2( A + A0 2 (A + A 2 ) 
x exp[ — (Ai + A 2 + 2A)]- 



<i ' - ' /J (l + T A 1 )V2(i +ToA 2 )V2 
Ax A 2 2A .o 



_^ 9T 2 r Ai(l + AQ A 2 (l + A 2 ) 2A(1-A) 

+ZI ° ^I + TqAO^ (l + ToA 2 ) 2 + (1-T A) 2J } ' W 



The four integrals must be done numerically. Eq. fl24f) coincides with the result obtained by 
straightforward application of the supersymmetry approach. We observe that the integrand 

9 



in Eq. (|24|) is semi-positive definite. Therefore, L m decreases monotonically with increasing 



Tiv/c?, tending to zero as Y^/d — > oo. This is what we expect on physical grounds JT2]. For 
Tn /d = 0, we have I in = (Y^/Ys)/(l + Y^/Ys). This follows from I m = 1 and from Eqs. (f| 
0) and (|H). 

The appearance of the exponential factor in Eq. (^4|) can be visualized in yet another way. 
Instead of the replacement e — > — iY^, we could have put e = but kept in Eq. (8.10) of 



Ref. [13| the product over transmission coefficients 



n (1 ~ T * A) (25) 

V(l + TAi) 1 / 2 (l + TA 2 ) 1 /2 • l ) 

We could have argued that this product accounts for both, coupling to the SD band via T , 
and for El decay into a large number of open decay channels described by the transmission 
coefficients 7] with I ^ 0. Owing to the weakness of the electromagnetic force, we would have 
T[ <^ 1 for all I ^ 0, although the sum X^o^i ma y be significant. Excepting / = 0, we could 
then approximate the product (|25| ) by exp(— (l/2)(J]/7])(Ai + A2 + 2A)). The prime indicates 
that the term with / = is omitted. Comparison of this exponential with the one in Eq. ( p^]) 
shows that Yl[Ti = 2irYN/d. This is a very satisfactory result. It is identical to the standard 
relation connecting decay width and sum over transmission coefficients in the theory of nuclear 
reactions, cf. Ref. |H|. Hence, we identify the total transmission coefficient T/v for El decay as 



T N = 27rT N /d . (26) 

We note that in nuclei with mass A ~ 190 where T^/d ~ 10~ 2 we have T/v ~ .05. This is a 
rather small value. We will see in the next Section that owing to this small value, decay of the 
ND states back into the SD channel is not altogether negligible, in contrast to the claim made 



in Ref. [[12 



In summary, we have described how to generate an analytical expression for I in . As a 
by-product, we have seen that this quantity depends on the input parameters T s , T N , 1^ and d 
of our model only via T and T N , see Eqs. (^3]) and (^). Integration over energy E will reduce 
the dependence on the dimensionless parameter To to a dependence on the only surviving 
dimensionless combination r^/r^ of parameters contained in To. Hence, I in is a function of 
the two dimensionless parameters I^/T^ and Y^/d. This conclusion applies not only to I in 
but likewise to all higher moments of A Q0 (E). Indeed, the supersymmetry approach |L3| yields 



formally identical results for these higher moments, except for the fact that the dimension of the 
supermatrices is increased. Thus, we have shown that the entire probability distribution of 7i n 
depends on the input parameters of our model only via the two dimensionless variables I^/Ts 
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and Y^/d. This is a non-trivial result for two reasons. First, from the four input parameters 
V s-, Ttv, r^ and d of the model, we can construct three independent dimensionless variables. 
The probability distribution of I iri depends on only two of them. Second, both dimensionless 
variables are given by the ratio of an electromagnetic and a nuclear quantity. Hence, both 
depend on the ratio of the fine structure constant and the strength of the strong interaction. 

5 Results 

We have calculated the fourfold integral in Eq. (|24]) numerically. We add the term L in , cf. 
Eq. (|18D , and denote the result by (/). We have also simulated the model of Eqs. ([I], 0, (|) 
numerically. This was done by drawing the matrix elements Vqj from a Gaussian distribution 
centered at zero with variance v 2 . The energies Ej were taken from an unfolded GOE spectrum 
with E located in the center of the semicircle. Typically, we used matrices of dimension 
K = 100 or bigger, and we calculated N = 10 4 or more realizations. The calculations were 
simplified by using for Aqq(E) the expression 

Ax)(£) = S —£ vv ■ (27) 

E-E 0+ iT s /2-^ l E J^ N/2 

We note that in the simulation, we calculate the total intraband intensity I in without introduc- 
ing the distinction between I m and I m . We used (I) to test the results of the simulation. 
The width of the probability distribution of I in was estimated as follows. With I(n) the value 
of /jn obtained in the n th realization (n — 1, . . . ,N), two sets labelled s$ with i — 1,2 were 
formed depending on whether I(n) < (I) or I(n) > (I), respectively, each set containing iVj 
realizations labelled \ii = 1, . . . , iVj. For % = 1, 2 we have calculated 



O; 



I Ni 






The results are shown in Figs. 2(a) to 2(f). Cases (a) to (f) correspond to six different 
choices of the ratio T^/d as indicated in the Figures. In all cases shown, the abscissa gives the 
ratio T^/Tg on a logarithmic scale. The top panels show I in as calculated from Eq. (|2~3]). We 
note that J in decreases monotonically with increasing Y^/d. At the same time, the maximum 
of Jjn shifts towards smaller values of r^/Ts. We note, moreover, that for Y^/d = 10~ 2 , the 
peak value is ~ 0.6 and, thus, not small compared to unity. The decay of the ND states back into 
the SD band is negligible only for overlapping ND states, i.e., for Y^/d > 1. The bottom panels 
show as solid lines the average values I in of the total intraband intensity obtained by adding I, 



■av 

m 
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(Eq. flT8|)) and I m (top panels). All curves decrease monotonically with increasing T^/Ts- 
This is physically plausible. The value of r^/r^ where I[ n equals 1/2 shrinks from »5x 10 4 
for T^/d = 10~ 6 to ~ 1 for T^/d = 10. The bottom panels also show as error bars the widths 
of the distributions as estimated by the quantities <7j in Eq. (|28|) . Not surprisingly the errors 
are biggest when the contribution of I m to I in is largest and shrink with increasing size of 
T^/d. The error bars reflect directly the role played by individual ND states in the El decay: 
For T^/d ^> 1, the ND states overlap strongly, their positions are irrelevant, and the errors 
are small. The opposite situation prevails for T^/d <C 1 where the location of the ND states 
closest to the SD state |0) is of tantamount importance. 

In order to make our work useful for the analysis of experimental data, we now present 
two fit formulas which approximately reproduce the relevant behavior of the average intraband 
decay intensity / in . We fit the curves for I in shown in the top panels of Figs. 2(a) to 2(f) and 
find 

—flue , ,r* Nn9mN r (0.4343 in(rVr 5 ) - 0.45(^)- ai303 )2 , N 

/ in fluc = (1 - 0.9139(-^f 2172 ) x exp{-^ (r//V a1477 } ' ^ 

We emphasize that this formula is not based on any theoretical arguments and presents the 
result of an approach based upon trial and error. In Fig. 3, we show a comparison between the 
fit formula fl29"|) (dotted lines) and the calculated values for I in (solid lines). We observe that 
the decreasing parts of the curves for I[ n (which are particularly relevant for the experimental 
determination of ]N-) are particularly well reproduced. For fixed values of T^/d, the difference 
between the fit value and the calculated value is, for any value of T^/r^, never bigger than .07 
and lies well within the limits of uncertainty defined by the error bars in the lower panels of 
Fig. 2. 

For the lower panels of Figs. 2(a) to 2(f) we fit the values of r^/r^ for which, at a given 
value of T N /d, the average intraband decay intensity I in assumes the value 1/2. We obtain 

rVr 5 = i + o.5429(r^/t/)-°- 8947 . (30) 

Formula fl3"0|) is in good agreement with the calculated values of T^/T s , see Fig. 4. 

6 Comparison with the Approach by Vigezzi et al. 

The experimental results for decay out of a SD band have frequently been compared with the 
approach developed by Vigezzi et al. Therefore, we investigate in this Section how the approach 
and the results of Ref. [f| compare wih the theory developed in previous Sections of this paper. 
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(c) r N /d=io" 
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(e) r N /d=1.0 



(f) r N /d=io.o 
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Figure 2: The ensemble average I m (top, abbreviated as If) and the total average intraband 
intensity I m with error bars due to the statistical uncertainty (bottom, abbreviated as /) versus 
logarithm (basis ten) of the ratio r^/T^ for several values of T^/d. 
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Figure 3: Comparison of the fit formula Eq. ( p9l) (dotted lines) with the ensemble-averaged 
intraband decay intensity (solid lines) J in (abbreviated as Ip) for several values of Y N /d. 
The basis of logarithm is ten. 
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Figure 4: Comparison of the fit formula Eq. (|30D (solid line) with the calculated values of 
T^/Ts (circles) for which the average intraband decay intensity attains the value 1/2. The 
basis of logarithm is ten. 
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The model which serves as starting point of the approach of Ref. || is identical to the one 
used above, see Eqs. (|I|, 0, [3], |J). The formula actually used by Vigezzi et al. to calculate 
the probability J out for decay out of the SD band is, however, not really derived from this 
model. It is rather based on physically plausible and intuitive reasoning which we summarize 
as follows. Let \m) with m = 1, . . . , K + 1 denote the (K + 1) eigenstates of the Hermitian 
matrix H defined in Eq. ([!]). We note that this matrix does not contain the decay widths 
Y$ and Tjv which appear in the width matrix, Eq. (fj), and in the effective Hamiltonian H c g, 
Eq. (Q). Let c m = (0\m) be the amplitude with which the SD state |0) is admixed to the state 
\m). It is argued that the width Ys{m) for decay of the state \m) into the SD band is given by 
Ys{m) = |c m | 2 r£. The width Y^f(m) for El decay of the state \m) is correspondingly written 
as Yn(i7i) = (1 — |c m | 2 )rAr. Finally, it is assumed that E2 decay out of the next-higher state 
in the SD band populates the state \m) with probability \c m \ 2 . The intensity I out for decay out 
of the SD band is then given by summing over all states \m) as 

pig _ V^ I 12 U " \ C m\ JMV /oi\ 

V (1 - \c m \ z )Y N + |c m |"T s 

We have added a superscript vig to identify the origin of this formula. We note that according 
to Eq. (|3lD, the probability distribution of J^f and, therefore, also that of the intraband E2 
intensity I^ g = l — I^, depend on two dimensionless variables: The ratio T^/Tat which appears 
explicitly in Eq. (pT|), and the ratio Fyd which determines the statistical behavior of the mixing 
parameters \c m \ 2 . We note that this parametric dependence of the Vigezzi model differs from 
the one characterizing the exact theory of Section [| where the relevant parameters are T^/d 
and T^/Ts- The reasoning behind Eq. ( pi]) leads us to expect that this equation renders a useful 
approximation to the exact result whenever T$ and Tn are sufficiently small (case of isolated 
resonances, T s and T N <C d) . This is also suggested by the observation that J^f is independent 
of the value of the fine-structure constant, a result which is not physically plausible. The 
worrisome aspect is that analysis of the data using the Vigezzi approach yields values of T^ 
which are about two orders of magnitude smaller than T^ 0, H, |9], [10|], thereby putting the 
entire approach into question. This is in fact what prompted the work of Ref. [TI| as well as the 
present investigation. By comparing numerical results obtained for 7j^ lg with those generated 
from the model of Eqs. fll], 0, |3|, f|), we now display the limits of validity of the approach of 
Vigezzi et al. 

Fig. 5 shows the average total intraband decay intensity calculated from the theory of 
Section [| (solid lines) and the result of the approximation ( |5T] ) (dotted lines) versus T^/Ts for 



six choices of Y^/d. We see that the difference between the approximation by Vigezzi et al. 
and the full theory becomes significant only at values of Y^/d w 1 which lie well outside the 



range which is relevant for the data of Refs. [0, ||, || [10| as well as the expected range of validity 
of the Vigezzi approach. The plots in Fig. 6 use the parameter Fyd of the Vigezzi model to 
define the abscissa and show similar behavior. Further insight is provided by considering the 
case 1 ^> r-L/Ts ^> T^/d which combines weak mixing between SD and ND states with strong 
fluctuations. In this case, a modified perturbation treatment yields to lowest order 



'7r rJ-r^r 



out V 2 T s d 



for i > rVr s > r N /d . (32) 



A related formula was given by Vigezzi et al. [|j . The predictions of this formula are shown as 
dash-dotted lines in Fig. 6. We note that whenever the condition of validity 1 ~^> T^/r^ ^> T^/d 
is met, the predictions of Eq. ( |32|) agree very well with the exact result. 

Can we determine the limits of validity of the Vigezzi approach also from theoretical argu- 
ments? As noted above, inspection of Eq. (|3~T|) suggests that the approach is correct to lowest 
non-vanishing order in both Y$ and T^ (isolated resonances). Moreover, it seems to account for 
all orders of Vqj. The restriction to non-overlapping resonances becomes obvious by considering 



-flue 



the case T N /d > 1. Here, I in is negligible, and 7 out is simply given by (r^/r s )/(l + (r^/T^)). 
This result is obviously and not surprisingly inaccessible to Eq. (|3TD. For the case of strong 
coupling T^ 3> T s , it simplifies to J out ~ 1 — (Ts/T^-). This is in contrast to the expression 
(1 — dT $ / (nF^F n)) derived in Ref. |4j] for the same regime. Aside from this rather trivial point, 
a better understanding of the domain of validity of Eq. ( |3T| ) is obtained by considering two 
limiting cases, the limit r^ — > with T^ ^ fixed (limit A), and the limit T^ —>■ with T^ ^ 
fixed (limit B). For limit A, Eq. (p4|) shows that J? uc vanishes as (r^) 2 . Keeping only linear 
terms in T^, we have from Eqs. (|17|) and ([18]) that I in = (1 + r^/r^) -1 ~ 1 — T^/Tg and, by 
Eq. (§), that 7^ « T^Ts- This agrees with Eq. fllj). For limit B, we recall that for T N = 0, 
the right-hand-side of Eq. ( P^ ) must equal (T^/Ts)/(1 + (T^/Ts))- (Admittedly, this is not 
immediately obvious from Eq. (|24]) ). The lowest-order correction to this result is obtained by 
expanding the exponential in Eq. (|4]) in powers of T^/d. The first-order term is negative and 



yields for 7 out a positive term linear in T^/d. While Eq. fl3T|) is in keeping with limit B, this 
equation is not consistent with limit A: I^f never becomes independent of T^, no matter how 
small VK This confirms the belief stated in the Introduction that the approach of Vigezzi et al. 
applies only for sufficiently large values of I^. To quantify this statement, we have calculated 
the terms of next order in the perturbation expansion of Section |3|. We find that, for Tn -C d, 
these are of order (T^/Ts) 2 {d/TN). Limit A is excluded if the second-order terms are at least 
of the same order as the first-order ones, i.e., whenever T^/Ts > T 'n/cL or T^ < d^T^/Ts). Vi- 
olation of this condition accounts for the deviations between the exact results and the Vigezzi 
approach displayed in Fig. 6. 
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Figure 5: Comparison of the ensemble-averaged total intraband decay intensity (solid lines) 
Jin (abbreviated as /) with the prediction of the Vigezzi approach (dotted lines). The basis of 
logarithm is ten. 

20 



•3-2 



/ 
/ 
/ 
/ 
/ 
/ 




s*' 


r s /r N =o.oi,d=i.o 

Vigezzi model 

full theory, T N = 0.001 
full theory, T N = 0.1 

— - full theory, T N = 10.0 
Formula (32) 



^ -1 



^50-2 
o 



,.' 




~~-~~~^' — iJi!SW 


*■' ''^sl**** 


/• 


, ^fi^ / 




-■ '^/ ^ 




/' 






/ 


r s /r N =2.o, d=i.o 

Vigezzi model 


/ 


full theory, Tn= 0.01 


/' 


full theory, T N = 0.1 




— - full theory, T N = 1.0 


/ 


Formula (32) 





^. <^/ 






z***^ 


"•-2 






'Ho 
2-4 


/ 


r s /r N =ioo.o, d=i.o 


/ 


Vigezzi model 

full theory, F N = 0.00001 




/ / 


full theory, F N = 0.001 




/ 


— - full theory, l' N = 0.1 


-6 


/ 


— — Formula (32) 



iog(ri / d) 



Figure 6: The logarithm of the ensemble-averaged decay intensity I out out of the SD band 
versus the logarithm (basis ten) of the ratio T^/d for three values of r^/TV- We show predictions 
of the Vigezzi approach (solid lines), results of the full theory for several values of T^, and the 
prediction of Eq. fl32|). 
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In conclusion, we see that the approach by Vigezzi et al. is subject to two constraints. The 
obvious one is that it deals with isolated resonances. This implies T^ <C d. The second, less 
obvious one is due to the constraint Tjy < d(T^/Ts)- 

7 Summary 

In the present paper we have calculated the ensemble average and properties of the distribu- 
tion function of the intraband E2 decay intensity for a statistical process leading to decay out 
of a SD band. We have shown that the entire distribution function depends only on the two 
dimensionless ratios T^/Ts and T^/d. Writing the intraband intensity as the sum of two terms, 
given in terms of the average decay amplitude and of its fluctuating part, respectively, we have 
shown that / dominates for large values of T^/d > 1, while both the average fluctuating part 
/ " and the fluctuations around it attain maximum values for small values of T^/d. We have 
proposed two fit formulas. One permits an estimate of the average intraband decay intensity 
and the other, an estimate of the value of T^/Ts for which the average intraband decay intensity 
attains the value 1/2. We have compared our results with those of Vigezzi et al. We have shown 
that the latter approach is an approximation to the exact theory, and we have established its 
limits of validity. For practical purposes, the approach offers a useful approximation whenever 
V N /d< 1CT 1 . 
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